Radial mode dependence of optical beam shifts 
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It is known that orbital angular momentum (OAM) couples the Goos-Hanchen and Imbert-Fedorov shifts. 
Here, we present the first study of these shifts when the OAM-endowed LGi p beams have higher-order radial 
mode index (p>0). We show theoretically and experimentally that the angular shifts are enhanced by p while 
the positional shifts are not. © 2012 Optical Society of America 
OCIS codes: 240.3695, 260.5430. 



A bounded beam of light may experience a displacement 
and/or a deflection with respect to the geometric optics 
prediction upon reflection [1,2]. These shifts are known 
as the Goos-Hanchen (GH) [3] and the Imbert-Fedorov 
(IF) [4,5] shifts. A shift in the plane of incidence indicates 
that the beam underwent a GH shift while a shift that 
is perpendicular to the plane of incidence connotes an 
IF shift. Moreover, these shifts are further distinguished 
because of their positional and angular parts [1]. The 
positional shift is the displacement of the center point of 
reflection while an angular shift happens when the beam 
deflects and propagates with a slight deviation from the 
law of reflection. The angular shift is responsible for the 
increasing beam excursion from the geometric optics pre- 
diction with propagation. All these shifts are in the order 
0.1-10 wavelengths for the positional shifts and micro- to 
milliradians for the angular shifts. 

These shifts may occur simultaneously or separately. 
Their presence depends on the polarization of the inci- 
dent beam [1-5], the index gradient seen by the beam 
[6,8-10], the beam's divergence (opening angle of the 
beam) [1,7,11] and, as emphasized recently, on the modal 
structure of the beam [12-16]. In works that involve 
modal structures, the Laguerre-Gaussian beam is used 
as the definitive beam because of its orbital angular mo- 
mentum (OAM). The OAM of light couples the four 
beam shifts (the positional GH and IF shifts and their 
angular counterparts). First described in [12], we subse- 
quently proposed a matrix which was validated with ex- 
periments, that neatly summarizes this mixing [13]. But 
since we considered only the case of a LG radial mode 
index that is zero (p = 0), that matrix is only explicitly 
dependent on the azimuthal mode index, I. The most 
important result of [13] is the measurement of the OAM- 
induced positional GH and IF and the OAM-affected an- 
gular GH and IF shifts for £ = —1,1. We extended the 
measurement for \£\ > 1 in [17]. The shifts are found to 
be linearly enhanced by the OAM. 

In this Letter, we answer the question to what degree 
the complete spatial mode of a beam affects these shifts. 



We first do theoretical calculations and extend our previ- 
ous matrix formulation to the case of p > . Next we re- 
port experiments on shifts associated with higher-order 
mode Laguerre-Gaussian modes of different azimuthal 
and radial mode indices. 

Laguerre-Gaussian modes are a set of solutions of the 
paraxial wave equation of the form given by 
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where x = r cos cf>, 
tal beam waist, L 



y = r sin <fi, wq is the fundamen- 
= fcoiUo/2 is the Rayleigh length, 



w (z) = u>oy 1 + i z /L) is the beam waist at z, kg is 
the center wave vector, and L^(x) denotes the general- 
ized Laguerre polynomial of order (£, p) with p > and 
I as integers. This beam is characterized by the indices 
\l\ and p , where £ is the azimuthal mode index and 
p+ 1 is the number of radial nodes [18]. There has been 
little interest in the index p ; mainly for use in optical 
manipulation and trapping [19-22]. However, higher LG 
modes (p,\£\ > 0) are increasingly being noticed, for in- 
stance in gravitational wave detection. As a matter of 
fact, some tabletop experiments are being done to test 
their suitability [23-25]. 

The theoretical derivation of the effect of the radial 
mode index follows from [1, 13]. The incident and re- 
flected beams are assumed to be Gaussian. The incident 
beam is then decomposed into plane wave components 
and the Fresnel reflection coefficients are applied to the 
s and p polarization components of the wave, respec- 
tively. The shift is obtained by summing all the reflected 
plane wave components, and taking the centroid of the 
intensity distribution. 
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Following the notation in [13] and noting that we are 
using a median detector [26,27], we arrive at a matrix 
which is explicitly dependent on both £ and p , given by, 
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where Aq H and A IF are the dimensionless positional 
GH and IF shifts, respectively, and Qq H and Q° IF are 
the dimensionless angular GH and IF shifts, respectively, 
for a TEMqq beam. These are given by, 
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Figure 1. (color online) The experimental set-up used 
to measure the polarization-differential beam shifts as 
a function of the incident angle. LCVR is the variable 
retarder. SWG is the square wave generator. QD is the 
quadrant detector. LI, L2, and L3 are lenses. More de- 
tails are in the text. 
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with w s/p = R 2 s/p a 2 s/p / (R 2 p al + R 2 s aj), and r s/p = 
R s / P exp (itps/p) the Fresnel reflection coefficient evalu- 
ated at the incident angle 9i and a s / p the electric field 
components. The shifts x and y are the sum of two con- 



tributions and are given by kgx = A 



GH 



- (z/ L)Q GH 



and 



k y = A IF + (z/L)Q IF , respectively. 

Finally, the matrix element £ (£,p) calculated for a me- 
dian detector is £ (£, p) = 21+ / /_ where 
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It follows from (2) that the positional shift is not af- 
fected by the p index. Only the angular shift components 
are influenced by it. Moreover, the matrix elements M22 
and M44 are the same. This means that the enhance- 
ments of the angular GH and IF shifts are the same. 

Our experimental setup is shown in Fig. 1. We use 
an open cavity HeNe laser (A = 632. 8nm) and placed 
two crossed copper wires (d ~ 90/im) inside the cav- 
ity. The wires force the cavity to emit a higher-order 
Hermite-Gaussian mode (HG m ,n)- The integer mode in- 
dices m and n can be dialed by properly shifting the 



intracavity wires. With an astigmatic mode converter, 
the HG m . n modes are then converted into LGi tP modes 
with transformations £ = m — n and p = min(m,n) 
as described in [28,29]. LI is the mode matching lens 
while lens L2 collimates the beam. The beam then passes 
through an LC variable retarder (LVCR) which sets the 
beam's polarization state. The LCVR (Meadowlark) is 
driven by a square wave generator at 2.5 Hz. This wave 
generator also provides the reference wave of the lock-in 
amplifier (EG&G 5210). The positional shifts are meas- 
ured with a collimated beam while the angular shifts are 
measured with a beam focused by lens L3 (f=70 mm). 
The displacement due to the angular shifts increases 
with the opening angle of the beam 9q, as Oq/2 where 
f5 = 2/(fc wo) [11,30]. 

The polarization differential shifts are measured be- 
tween 45°/ — 45° polarization states. These states are 
specifically chosen because: 1) the dimensionless angular 
IF shift Q IF that is responsible for the Aq H and Q e IF 
is maximum at these states; and 2) the contribution of 
A°j ff , Aj F , and Q% H either have values equivalent to 
zero or cancel to zero upon measurement of the differen- 
tial shift. These enable us to isolate the matrix elements 
M12 and M22 in (2). 

After the polarization conditioning, the beam is made 
to reflect at different incident angles 6i on a BK7 prism 
(n = 1.51) mounted on a 9 — 29 rotation stage. The 
position of the reflected beam is measured by a quadrant 
detector (QD) whose signal is fed to the lock-in amplifier. 
These are done for different combinations of £ and p. 

Figure 2 shows the measured matrix element M12 as 
function of £ and with different values of p . The meas- 
ured values collapse to the Mi 2 = — 2£ line. This is 
correctly predicted by our theory. The matrix element 
changes linearly on £ but is independent on p. We have 
experimentally verified that the positional GH shift is 
not affected by the radial mode index (not shown). 

The measured M22 is shown in Fig. 3. The theoretical 
predictions are calculated from £ (£,p) for different £ and 
p values. For p > we obtained greater enhancements. 
This confirms our theoretical prediction of the effect of 
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Figure 2. (color online) Matrix element for different 
values of I and p. The solid line is Mn = —21. 
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Figure 3. (color online) Matrix element M22 for different 
values of £ and p. The solid line is £ (£,p). A discussion 
of £ (£,p) is in the text. 

p on the angular shifts. At p = 0, the angular IF shift 
linearly increases as 1 +2|^|. Furthermore, had we used a 
centroid detector the increase of the angular shift would 
have been £ (£, p) = 2p+ \£\ + 1, a linear increase with p. 

In conclusion, we have theoretically described and ex- 
perimentally measured the positional GH shift and an- 
gular IF shifts for higher-order LGi iP to determine the 
role of mode indices I and p. Specifically, we have ob- 
served that only the azimuthal index I affects the po- 
sitional shift while the angular shifts are enhanced by 
both the radial p and the azimuthal £ indices. Higher- 
order LGi tP beams form a complete basis set for paraxial 
light beams contrary to LGi p= o [18]- Hence with proper 
mode decomposition, beam shifts associated with arbi- 
trary modes of light can be predicted based on our the- 
oretical and experimental results. 

This work is supported by the Foundation for Fun- 
damental Research of Matter (FOM) and the European 
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Program 255914 Phorbitech. 
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